Transport Anisotropy as a Probe of the Interstitial Vortex State in Superconductors 

with Artificial Pinning Arrays 
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We show using simulations that when interstitial vortices are present in superconductors with 
periodic pinning arrays, the transport in two perpendicular directions can be anisotropic. The degree 
of the anisotropy varies as a function of field due to the fact that the interstitial vortex lattice has 
distinct orderings at different matching fields. The anisotropy is most pronounced at the matching 
fields but persists at incommensurate fields, and it is most prominent for triangular, honeycomb, 
and kagome pinning arrays. Square pinning arrays can also show anisotropic transport at certain 
fields in spite of the fact that the perpendicular directions of the square pinning array are identical. 
We show that the anisotropy results from distinct vortex dynamical states and that although the 
critical depinning force may be lower in one direction, the vortex velocity above depinning may also 
be lower in the same direction for ranges of external drives where both directions are depinned. For 
honeycomb and kagome pinning arrays, the anisotropy can show multiple reversals as a function 
of field. We argue that when the pinning sites can be multiply occupied such that no interstitial 
vortices are present, the anisotropy is strongly reduced or absent. 

PACS numbers: 74.25.Qt 



I. INTRODUCTION 

Vortices in superconductors interacting with artificial 
arrays of periodic pinning exhibit a wide range of coni- 
mensurability and dynamical effects that can be observed 
readily in critical current, transport, and other bulk mea- 
surements P, U, S i, S, i, 0, H, la, [i3 • Advances in lithog- 
raphy techniques permit the creation of pinning arrays in 
which the size, shape, and composition of the individual 
pinning sites and the glob al g eometry can be well con- 
trolled 0, 1 i, i, H, 0, a i, [Tfl [m [H, [H, [il . Commen- 
surability effects in these systems occur when the num- 
ber of vortices equals an integer multiple of the number 
of pinning sites, resulting in peaks or anomalies in bulk 
measurements as a function of field. At the first matching 
field, there is one vortex per pinning site, and as the field 
is further increased, additional vortices can be located 
either at the pinning sites in the form of pinned multi- 
quanta vortices or in the inter- 
stitial regions between the pinning sites. The interstitial 
vortices can be effectively pinned by the repulsive interac- 
tions from the vortices at the pinning sites, which create 
a caging potential [1 i, HI, Ho, [HI, [22, [H IH] . It is also 
possible for mixed vortex pinning to occur in which the 
first few matching fields have only pinned multi-quanta 
vortices until the pinning sites are saturated, while for 
higher matching fields the additional vortices are located 
in the interstitial regions ^, ^, JJ., 15,, J^, JJi] . Conversely, 
it is also possible that interstitial vortices appear at the 
lower matching fields, but that as the vortex- vortex inter- 
actions increase at higher matching fields, multi-quanta 
vortices will begin to form at the pinning sites jl9| . 

Interstitial vortex lattice crystals in square periodic 
pinning arrays have been observed directly with Lorentz 
microscopy, which revealed that there are several dis- 
tinct types of interstitial vortex structures that have 



symmetries different from that of the triangular vortex 
lattice The same types of vortex structures have 
been produced in simulations of square pinning arrays 
(20I [21I I , while simulations have also shown that similar 
vortex structures can form in triangular [20], rectangu- 
lar 24 [2^ , honeycomb [26| , and kagome pinning arrays 
I26I 27 1 . Other numerical works indicated that a rich va- 



riety of composite lattices with m ultip le and interstitial 
vortex configurations are possible Oll^l and that new 
types of interstitial vortex configurations can occur for 
arrays of antipinning sites [2^ . 

Vortex imaging experiments provide direct evidence 
for both multi-quanta vortex pinning and the formation 
of ordered interstitial vortex lattice structures il6| . 
Anomalies at matching fields found in bulk measure- 
ments occur for both multi-quanta vortex pinning and 
interstitial vortex pinning, so without direct imaging it 
can be difficult to determine whether multi-quanta or in- 
terstitial vortex pinning is occurring [isj . In some cases, 
the presence of interstitial vortices can only be inferred 
from the shapes and characteristics of the current volt- 
age curves or from phase locking experiments [2^ . It 
would be highly desirable to identify additional clear sig- 
natures in transport measurements that can distinguish 
between interstitial vortex pinning and multi-quanta vor- 
tex pinning and that can also reveal the types of vortex 
lattice symmetries that are present. 

An anisotropic response was recently measured for the 
critical current applied in two perpendicular directions 
to a triangular pinning array in recent experiments and 
simulations (sol . [3]| . The experiments were performed 
on several different samples and the same anisotropic re- 
sponse appeared in each one, while the anisotropy ob- 
served in the simulations agreed with that seen in the 
experiments, suggesting that the behavior is due to dis- 
tinct intrinsic features of the vortex dynamics. The sim- 
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ulations show that the vortex flow patterns are different 
for the two directions of applied current, which could 
account for the anisotropic response. The anisotropy is 
particularly pronounced at the second matching field but 
is absent at the first matching field, which suggests that 
depinning of the interstitial vortices is responsible for 
the anisotropy. Interestingly, for lower temperatures the 
experiments showed that the critical current anisotropy 
vanished at both the first and second matching fields 
but persisted at the third matching field, suggesting that 
multi-quanta vortex pinning occurs at the second match- 
ing field at low temperatures. These results indicate that 
the presence of anisotropy can be a useful way to probe 
the interstitial vortex state and the type of vortex order- 
ing that occurs at different matching fields. 

The experiments and simulations of Refs. [33,|3l| only 
examined a triangular pinning array up to the third 
matching field. In this work we study anisotropic trans- 
port for a much wider range of fields and system pa- 
rameters for triangular, square, honeycomb, and kagome 
pinning arrays. For the triangular array we find a vari- 
ety of distinct anisotropic behaviors due to the differing 
symmetries of the interstitial vortex lattice at different 
matching fields. For example, at certain matching fields 
the vortex lattice is disordered and the anisotropy van- 
ishes. For honeycomb and kagome pinning arrays, an 
even richer anisotropic behavior occurs due to the for- 
mation of vortex molecular crystal states [1^ with addi- 
tional rotational degrees of freedom, resulting in a series 
of reversals in the anisotropy as a function of field. Re- 
markably, we find that it is also possible for square pin- 
ning arrays to show anisotropic transport even though 
the two perpendicular directions of the pinning lattice are 
identical. This occurs at certain matching fields where a 
triangular vortex structure forms, such that one driv- 
ing direction is oriented with the easy-shear direction 
of the vortex lattice. Anisotropic transport measure- 
ments for vortices in periodic pinning arrays have already 
been been shown to be experimentally feasible, and ex- 
periments have been performed in which the current is 
injected in two directions for samples with rectangular 
pinning arrays or asymmetric pinnin g sh apes, revealing 
anisotropic depinning thresholds [3^. l33l. [33| . Although 
our work is focused on superconducting vortices, our re- 
sults should be general to the class of systems of particles 
interacting with periodic substrates where both pinned 
and interstitial particles are present. Examples of this 
type of system include colloidal particles in periodic pin- 
ning [1^, |3^ and charged metallic balls [s^l ■ 



II. SIMULATION AND SYSTEM 

We consider a two-dimensional system with periodic 
boundary conditions in the x and y-directions of size 
L X L. The magnetic field B is applied out of the plane 
in the z-direction, while Ny vortices and Np pinning 
sites are placed within the system for a vortex density of 
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FIG. 1: The pinning site locations (open circles) for (a) a 
triangular pinning array, (b) a square pinning array, (c) a 
honeycomb pinning array, and (d) a kagome pinning array. 
The driving force is applied along y to determine the depin- 
ning force and along x to determine F^. 



riy — Ny/ LP' and a pinning density of — N^jl? . The 
matching field is defined as the field at which the 
number of vortices equals the number of pinning sites, 
= Ap. In Fig. [1] we show representative examples of 
the triangular, square, honeycomb, and kagome pinning 
geometries used in this work, and indicate the x and y 
directions along which current is applied. The initial vor- 
tex positions are prepared by simulated annealing with 
no applied drive, and then the vortex velocities are mea- 
sured in the presence of a driving force that is applied 
in the x direction. The driving force corresponds to the 
Lorentz force generated by an applied current, while the 
vortex velocities are proportional to the voltage response 
that would be measured experimentally. We repeat the 
simulation from the same initial vortex positions with the 
driving force applied in the y direction, and compare the 
velocity response and (Vy) and the critical depinning 
force FJ* and for the two driving directions. 

The time evolution of the vortex dynamics is governed 
by integrating A'^ coupled overdamped equations of mo- 
tion. The equation of motion for a single vortex i at 
position Ri is given by 



77- 
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(1) 



Here the damping constant r\ — 0Qd/27r^'^p7v, where d is 
the sample thickness, 77 is the coherence length, pn is the 
normal-state resistivity, and = h/2e is the elementary 
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flux quantum. The vortex-vortex interaction force is 

Fr=£/oifi(^)R.„ (2) 

wtiere Ki is tlie modified Bessel function, A is tlie Lon- 
don penetration depth, the unit offeree is /o = 0o/27rA'^, 
Rij = [Ri — Rjl, and Ry = (R, — Rj)/i?jj. The vortex- 
vortex interaction force falls off sufficiently rapidly that a 
cutoff can be placed at Rij = 6A. Use of a longer cutoff of 
Rij = 12A produces identical results. An additional short 
range cutoff is placed at Rij = 0.1 A to avoid a divergence 
in the force. The pinning sites are modeled as attrac- 
tive parabolic wells of radius Rp and strength Fp with 

= foFpR-^R^^eURp - R^^)/m^\ Here 
R^^-* is the location of pinning site k, Rik = jR^ — R[f''|, 
r(^) = (R,, - R[P^)/Rik, and 9 is the Heaviside step 
function. The pinning sites are arranged in a triangu- 
lar, square, honeycomb, or kagome array. The external 
drive Fjj — Fdx or Fj^ = Foy is a constant force that is 
uniformly applied to all of the vortices. The thermal 
force Ff is used during the simulated annealing pro- 
cedure and has the following properties: {Ff{t)) = 
and {F^{t)F[{t')) = 27jkBTSijS{t - t'). We decre- 
ment the temperature by 0.0002 every 1000 simulation 
time steps. After the initialization, the applied drive 
is imposed in increments of SFd =0.1 every 10'^ sim- 
ulation time steps. The velocity- force curves are ob- 
tained by averaging the velocity every 10^ simulation 
time steps: (Va) = N~^'^f''Vi ■ a, where a — x, y. 
Here Vj = dRi/dt. The critical depinning forces in the 
X and y directions, F^!^ and FJf, are determined by the 
criterion (Va) = 0.001. 

III. ANISOTROPY IN TRIANGULAR PINNING 
ARRAYS 

We first consider the depinning forces for the x and 
y-directions in the triangular pinning lattice illustrated 
in Fig. [Ha) with Fp = 0.85, Up = COSSS/A^, and 
Rp = 0.35A. In Fig.[2]we plot F^ and FJ' vs B/B^ along 
with a detail of the region from 1.0 < B/B^ < 5.0. In 
general F^ > F^ at most of the commensurate fields. 
The anisotropy at incommensurate fields is most pro- 
nounced for 1.0 < B/B^ < 3.0, as shown in the inset of 
Fig. [21 In previous simulations a similar anisotropy was 
observed for this range of fields [3l|. For B/B^/, > 3.0 
both F^ and F^ are small, so it is difficult to deter- 
mine whether the depinning is anisotropic at incommen- 
surate fields. At the matching fields B/B^ = 7, 9, and 
12, where the depinning forces are high, F|f and F^ can 
be measured accurately. In Fig. [3| we plot representa- 
tive velocity-force curves for driving in both the x and 
y directions at B/B^ = 2, 3, 4, and 5. There is a 
clear anisotropy in the depinning force with F^ > F^ at 
B/B^ = 2, 3, and 4 in Figs. [3);a,b,c), while at B/B^ = 5 




FIG. 2: Tlie depinning forces F^!' (liglit line) and FJ* (dark 
lines) vs B/B^ for a triangular pinning lattice with Fp — 0.85, 
Rp = 0.35A, and Up = 0.0833/A^ Inset: the same data 
showing a highlight of the region from 1.0 < B/B^ < 5.0. 
Open squares: F^; filled circles: . 
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FIG. 3: Velocity-force curves for driving in the s-direction, 
(Vx) (open squares), and in the y-direction, (Vy) (filled circles) 
for the triangular pinning lattice system in Fig. (2] (a) At 
B/B^ = 2, the slope dVy/dFo < dV^/dFo- (b) At B/B^, = 
3, dVy/dFo = dVx/dFD. (c) At B/B^ = 4, dV^/dFo < 
dVy/dFo, resulting in a crossing in the velocity- force curves 
near Fd = 0.1. (d) At B/B^ — 5, the depinning is isotropic. 



in Fig. [St^d), the depinning is isotropic. This can be 
seen more clearly in the plot of the anisotropy F^ / F^ at 
the matching fields, shown in Fig. (4] The anisotropy is 
largest for B/B^ = 2, 4, and 9, weaker for B/B^ = 3 and 
12, and essentially absent at B/B^ = 1, 5, 6, 7, 8, 10, and 
11. Corresponding to this, there are no peaks in the de- 
pinning force at B/B^ = 5, 6, 8, 10, and 11 in Fig. [2] At 
the fields with isotropic depinning F^ /F^ « 1, the overall 
vortex lattice is disordered, while at the other matching 
fields where anisotropic depinning occurs, ordered vor- 



FIG. 4: The ratio / vs B/B^ obtained from the system 
in Fig. 2. The dashed hne indicates F^ /F^ = 1. The depin- 
ning is strongly anisotropic at the matching fields B/B^ = 2, 
3, 4, 9, and 12. 



tex lattices form. Previous numerical work for triangu- 
lar pinning arrays at fields up to B/B,f, — 9 showed the 
same features in the critical depinning force as well as 
the existence of disordered lattices at the matching fields 
B/B^ — 5, 6, and 8 [13]. From a geometric construc- 
tion, it can be shown that a triangular vortex lattice can 
be placed on a triangular pinning lattice at the integer 
matching fields N — m? +12^ + nm, where n and m 
are integers [2^]. This predicts the formation of trian- 
gular vortex lattices at fields with N = 1, 3, 4, 7, 9, 
and 12, in agreement with our observation of a peak in 
Fc at each of these fields. The geometric construction 
does not predict the formation of a triangular lattice at 
B/Bfj, = 2; however, a strong matching peak appears at 
this field both in Fig. [2] and in the previous work [2^ . 
The peak at B/B^ = 2 occurs due to the formation of an 
ordered honeycomb vortex lattice structure, rather than 
a triangular vortex lattice. In general, we expect to find 
a peak in the critical current at fields where a triangular 
or other ordered vortex lattice structure forms. Figures [2] 
and [4] also show that although there is a peak in Fc at 
B/Bfj, — 7, where a triangular vortex lattice forms, there 
is no anisotropy for this field and F^ /F^ ~ 1. 

In order to explain the different degrees of anisotropy 
that appear at different fields, we analyze the vortex po- 
sitions and trajectories for driving in the x and y direc- 
tions. Figure [5l[a) shows the ordered honeycomb vortex 
lattice structure that forms at Fu = for B/B^ = 2 
on a triangular pinning lattice. When driven in the x- 
direction, the vortices can easily form slightly undulat- 
ing channels of flow that pass between the filled pinning 
sites, as illustrated in Fig. O^b). For driving in the y- 
direction, the pinned vortices act as barriers that prevent 
the formation of simple fiow channels, and the intersti- 
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FIG. 5: Vortex positions (black dots), pinning site locations 
(open circles), and vortex trajectories (black lines) for the 
triangular pinning lattice system in Fig. 2. (a) The vortices 
form a honeycomb configuration at Fd = and B/B^ — 2. 
(b) The vortex trajectories just above depinning for driving in 
the a;-direction at B/B^ = 2. The vortices channel between 
the pinning rows, (c) Vortex trajectories just above depin- 
ning for driving in the y-direction at B/B^ = 2. The moving 
interstitial vortices wind around the occupied pinning sites, 
(d) The vortices form a triangular configuration at Fd ~ 
and B/B^ = 3. (e) The vortex trajectories just above de- 
pinning for driving in the a;-direction at B/B^ — 3 show an 
ordered interstitial fiow of vortices between the pinning sites, 
(f) Vortex trajectories just above depinning for driving in the 
y-direction at B/B^ — 3. Here the flow is disordered. 

tial vortices can only move by making significant excur- 
sions into the x-direction, as illustrated in Fig. [5^c) for a 
drive just above depinning. As a result, a larger external 
force is required to cause depinning in the y-direction, 
and F^ > F^ ■ The onset of motion in the y-direction 
is very sharp, as seen by the jump in the velocity-force 
curve in Fig.[3ja), and after a brief initial period of disor- 
dered motion, the vortex flow quickly organizes into the 
pattern shown in Fig. [5]Jc). We observe similar motions 
just above depinning in the x and y-directions at the 
incommensurate fields for 1.0 < B/B^ < 2.0; however, 
the presence of vacancies in the honeycomb vortex lattice 
causes certain rows of interstitial vortices to depin at a 
lower value of Fu than at the commensurate fields. The 
velocity-force curves at B/B^ — 2.0 shown in Fig. [2Ia) 
indicate that dVy/dFn < dVx/dFu, meaning that the 
vortex velocity has a lower slope as a function of increas- 
ing drive in the y-direction than in the x-direction, even 
though the same number of vortices are moving for either 
direction of drive. This difference is a result of the fact 
that the vortex motion for y-direction driving has much 
larger excursions transverse to the driving direction since 
the interstitial vortices must go out of their way to pass 
around the pinned vortices. 

At B / B^ = 3 on a triangular pinning lattice, the over- 
all vortex lattice ordering at Fd = is triangular, as 
shown in Fig.[5Jd). The vortices move in slightly winding 
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channels upon application of a drive in the x-direction, 
as seen in Fig. Efe); these structures are similar to the 
channels that form at B/B^ = 2 in Fig.[5l^b). For depin- 
ning in the y-direction at B / B^ = 3, the pinned vortices 
again create a barrier to the formation of simple channels 
of interstitial vortex flow; however, unlike the ordered 
and strongly winding channels that form at B/B^ = 2, 
for B/B^ — 3 the vortices move in a disordered fashion 
just above depinning, as illustrated in Fig. ^i). An or- 
dered flow state similar to that shown in Fig. does 
not occur for B/B^ = 3 until a much higher value of 
Fjj is applied. This result indicates that in addition to 
the anisotropy in / ^ there are also strongly differ- 
ent vortex velocity fluctuation characteristics for driving 
along the x and y-directions at B/B^ — 3. Narrow band 
noise signatures should arise from the synchronized vor- 
tex motion that occurs for driving along the x-direction, 
while for driving in the y-direction the velocities are more 
random and a broad band noise signature should appear. 
In mode-locking experiments, where an external ac drive 
is imposed along with an applied dc drive, Shapiro type 
steps would appear for the ordered motion along the x- 
direction, while Shapiro steps would be absent for driving 
along the y-direction. Shapiro steps could be induced for 
both driving directions at B/B^ — 2 since both vortex 
flow patterns show synchronized motion; however, some 
of the characteristics of the Shapiro steps might differ for 
the two directions since the meandering of the vortices is 
distinct in the x and y-directions. 

A triangular vortex lattice with a single row of intersti- 
tial vortices between each row of pinning sites forms for 
B/Bfj, = 4 at Fo — 0, as shown in Fig.[6^a). Just above 
the depinning transition for driving in the x-direction, 
a portion of the interstitial vortices depin into the one- 
dimensional channels illustrated in Fig. [Sib). Here, one- 
third of the interstitial vortices remain immobile in the 
interstitial regions between pinning sites; these immobile 
interstitial vortices depin at a higher value of Fjj that 
is outside the range of driving forces shown in Fig. OJc). 
For driving in the y-direction, all of the interstitial vor- 
tices depin simultaneously and flow in winding channels 
around the occupied pinning sites, as shown in Fig.[S]Jc). 
The velocity- force curves shown in Fig.IS^c) indicate that 
although py > F^ for B/B^ = 4, dV^/dFo < dVy/dFo- 
As a result, the two velocity-force curves cross near 
Fd — 0.1. The slope dV/dFo is steeper for the y- 
direction driving since all the interstitial vortices take 
part in the motion, whereas for the cc-direction driving, 
only 2/3 of the interstitial vortices are moving. 

At B/Bij, — 5 and Fd = the vortex lattice is 
disordered, as seen in Fig. [Gfd), and the depinning is 
isotropic. Figure 6(d,e) shows that the same type of dis- 
ordered vortex motion occurs for depinning in both the x 
and y-directions. In general, we observe disordered flow 
states in both driving directions at the other matching 
fields where a disordered vortex lattice forms, including 
B/B^ = 6, 8, 10, and 11. 

At B/Bcji — 9, the depinning thresholds and velocity- 
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FIG. 6: Vortex positions (black dots), pinning site locations 
(open circles), and vortex trajectories (black lines) for the 
triangular pinning lattice system in Fig. (2] (a) At Fd ~ and 
B/Bif, = 4 a triangular vortex lattice forms, (b) The vortex 
trajectories just above depinning for driving in the a;-direction 
at B/B^ = 4. A portion of the interstitial vortices move in 
one- dimensional channels between the rows of pinning sites, 
(c) The vortex trajectories just above depinning for driving 
in the y-direction at B/B^ = 4. The interstitial vortices flow 
in winding channels around the pinning sites, (d) At Fd ~ 
and B/B^ — 5, the vortex lattice is disordered, (e) The 
vortex trajectories just above depinning for driving in the x- 
direction at B/ B^ = 5. The vortex flow pattern is disordered, 
(f) Vortex trajectories just above depinning for driving in the 
y-direction at B/B^ = 5. The same type of disordered flow 
pattern seen for x-direction driving appears for y-direction 
driving. 

force curves are very similar to those found for B jB^ ~ 4 
since the triangular vortex lattices that form at these two 
fields have the same orientation. In Fig. [Tl^a), the vor- 
tex configuration at -Fd = for BjB^ — 9 contains two 
rows of interstitial vortices between adjacent pairs of pin- 
ning rows, whereas at BjB^ = 4, Fig. [Slja) shows that 
there is only one row of interstitial vortices between each 
pair of pinning rows. Just above the depinning tran- 
sition for x-direction driving at B/B^ = 9, Fig. [7{b) 
indicates that the two rows of interstitial vortices flow 
in one-dimensional channels between the pinning rows 
while two interstitial vortices remain trapped behind ev- 
ery pinning site so that 3/4 of the interstitial vortices 
are moving. This is similar to the x-direction depinning 
that occurs for B/B^ = 4, where a single row of inter- 
stitial vortices flows between each pair of pinning rows 
and a single interstitial vortex is trapped behind each 
pinning site. Just above the depinning transition for y- 
direction driving at B/B^ = 9, shown in Fig. \T[^c), all 
of the interstitial vortices are depinned and a combina- 
tion of ordered and disordered vortex flow occurs. The 
velocity-force curves for driving in the x and y directions 
show a similar crossing at B/B^ = 9 as that illustrated 
in Fig.EJc) for B/B^ = 4. 

For B/B^ = 12, the vortex configurations and depin- 
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FIG. 7: Vortex positions (black dots), pinning site locations 
(open circles), and vortex trajectories (black lines) for the 
triangular pinning lattice system in Fig. 2. (a) At Fd = 
and B/B^ = 9, a triangular vortex lattice forms, (b) The 
vortex trajectories just above depinning for driving in the 
x-direction at B/B^ = 9. Two rows of interstitial vortices 
move in one-dimensional paths between the rows of pinning 
sites while a portion of the interstitial vortices remain pinned, 
(c) The vortex trajectories just above depinning for driving 
in the y-direction at B/B^ = 9. A combination of ordered 
and disordered flow occurs, (d) At Fd ~ and B/B^ = 12, 
a triangular vortex lattice forms which is aligned with the y 
direction, (e) The vortex trajectories just above depinning 
for driving in the a;-direction at B/B^ — 12. Three rows of 
interstitial vortices move between the rows of pinning sites, 
(f) Vortex trajectories just above depinning for driving in the 
y-direction at B/B^ — 12 show the existence of disordered 
flow. 



ning dynamics are similar to those seen for B / = 3. In 
Fig. md) the Fo = vortex configuration at B/ B^ — 12 
consists of a triangular lattice that is aligned with the 
y-direction. Since there is now an interstitial column 
of vortices that could depin and flow between adjacent 
columns of pinning sites, while there are no straight rows 
of interstitial vortices, it might be expected that < F^ 
for this field. Instead, Fig. g] shows that Fv/F^^ « 1.13, 
so the depinning is still easier in the x-direction than 
in the y-direction; however, the depinning anisotropy is 
much smaller than that which appears at B/B^ = 3. 
Vortex motion just above depinning at B/Bif, — 12 for 
driving in the a;-direction occurs in the form of three or- 
dered winding rows passing between each pair of pinning 
rows, as shown in Fig.[7l^e). This is similar to the motion 
at B / B^ — 3 shown in Fig. EJ^e) where one winding row 
of interstitial vortices moves between the pinning rows. 
For depinning in the y-direction at B/Bff, = 12, Fig.[7]Jf) 
shows that the vortex trajectories are disordered in a 
manner similar to that found for y-direction depinning 
at B/B^ ~ 3, as seen in Fig.OJf). 

At B/Bcf, = 7, a peak in Fc occurs as shown in Fig. [21 
but there is almost no anisotropy in the depinning thresh- 
olds, as seen in Fig. 21 FigurejSI^a) illustrates that at this 
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FIG. 8: Vortex positions (black dots), pining site locations 
(open circles) , and vortex trajectories (black lines) for the tri- 
angular pinning lattice system in Fig. 2 at B/B^ — 7. (a) At 
Fn = 0, a triangular vortex lattice forms that is not aligned 
with the X or y axes but is tilted at an angle 6 « —78° to 
the X-axis. This ground state is two-fold degenerate since the 
vortex lattice could have been tilted at the opposite angle, 
9 f» -1-78°, to the a;-axis. (b) The vortex trajectories just 
above depinning for driving in the a;-direction. All of the in- 
terstitial vortices are moving, (c) The vortex trajectories just 
above depinning for driving in the y-direction. In this case a 
portion of the interstitial vortices remains pinned. 



field for = 0, a triangular vortex lattice forms which 
is not aligned with either the x or y directions, unlike the 
configurations found at B/Bij, = 3, 4, 9, and 12, but is at 
an angle 9 « —78° to the x-axis. The vortices fiow in or- 
dered patterns just above depinning in both the x and y 
directions, as shown in Fig.|5Ub,c). Since a portion of the 
interstitial vortices remain pinned just above depinning 
for driving in the y-direction, whereas all of the inter- 
stitial vortices are flowing for driving in the x-direction, 
dVy/dFo < dV^/dFo at B/B^ = 7. The absence of the 
anisotropy at this fleld is likely due to the fact that the 
main symmetry axis of the triangular vortex lattice is not 
aligned with either the x or y directions, as is the case at 
the other matching flelds which show anisotropy. 

An interesting feature that we observe at B/B^ = 7 
which does not occur at the other matching fields we have 
investigated is a switching dynamics that can be induced 
within the pinned phase. This is illustrated in Fig. [51 
where in the initial state, shown in Fig. Hl^a), the vortex 
lattice is tilted at 6 —78° to the x-axis. The ground 
state is twofold degenerate since it could also have been 
aligned at 9 ^ +78° to the x-axis. When a small driving 
force is apphed to the ground state shown in Fig. |51^a) 
along -1-45° to the structural rearrangement of 

the vortices occurs accompanied by a domain wall that 
passes through the system from right to left, as shown in 
Fig.|nib,c). The vortices move by approximately one lat- 
tice constant in the y-direction as the domain wall sweeps 
past, and in the final state the vortex lattice is tilted at 
6 w +78° to the x-axis. The system can be switched 
back to —78° if an external drive is applied along —45° 
to the X-axis. We expect that for higher matching fields 
beyond the fields that we consider here, a switching ef- 
fect will be present for matching fields containing two or 
more degenerate ground states where the vortex lattice 
could be arranged in several possible ways. The applica- 
tion of an external drive will lower the energy of one of 
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FIG. 9: Vortex positions (black dots), pinning site locations 
(open circles), and vortex trajectories (black lines) for the tri- 
angular pinning lattice system in Fig. 2 at B /B^ = 7. A weak 
external driving force is applied at an angle in the positive x-y 
plane that induces a structural transformation in the vortex 
lattice between the two ground states, —78° shown in panel 
(a) and +78° shown in panel (d). During the transforma- 
tion, illustrated in panels (b) and (c), the interstitial vortices 
shift in the positive j/-direction by approximately one lattice 
constant and a domain wall traverses the system. 



the orientations. 



„ 0.1 




FIG. 10: (a) Depinning threshold in the a;-direction, (open 
squares), and j/-direction, (filled circles), vs Fp for the tri- 
angular pinning lattice system in Fig. 2 at B/B^ = 2. (b) 
The corresponding anisotropy ratio F'jljFc vs F^. H: the 
honeycomb ordering illustrated in Fig. Illf a): PP: the par- 
tially pinned phase shown in Fig. Illf b) where a portion of 
the pinning sites are unoccupied and the depinning is plastic; 
DT: the distorted triangular phase seen in Fig. Illf c') where 
most pinning sites are unoccupied and the depinning is elas- 
tic. The three different phases are visible as features in the 
anisotropy. (c) The depinning thresholds F^ and F^ vs Fp 
for the same system at B/B^ = 3. (d) The corresponding 
anisotropy ratio F^ /F^ vs Fp. 
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A. Effect of Pinning Strength 

We next study the effect of varying the pinning 
strength Fp on the anisotropy at the different fields. In 
Fig. [TUf a) we plot and F^ versus Fp for a sample 
with a triangular pinning lattice at B/B^ = 2, and in 
Fig. llOf b) we show the corresponding anisotropy ratio 
F]! / F^ . Since the pinning induces the honeycomb vor- 
tex lattice structure shown in Fig. [TIT a) at this field, as 
Fp decreases the vortices shift into configuration that is 
closer to a triangular lattice. To accommodate this shift, 
some of the pinning sites are vacated with decreasing 
Fp. For Fp > 0.85, the honeycomb lattice (H) structure 
is stabilized, all of the pinning sites are occupied, the 
anisotropy is fixed near F^ /F^ — 2.5, and the depinning 
thresholds and F^ do not change significantly with 
Fp. For 0.175 < Fp < 0.85, the pinning is not strong 
enough to stabilize the honeycomb lattice and a partially 
pinned (PP) lattice forms in which only some of the pin- 
ning sites are occupied, as illustrated in Fig. [TlTb) for 
Fp = 0.5. The depinning in this regime is still plastic 



FIG. 11: The vortex configurations for the triangular pinning 
lattice system at B/B^ — 2 from Fig. 10(a,b) at Fd ~ 0. (a) 
The pinned honeycomb lattice (H) at Fp — 1.0. (b) The 
partially pinned lattice (PP) at Fp — 0.5. A portion of the 
pinning sites are unoccupied and the vortices depin plastically, 
(c) The weak pinning regime at Fp = 0.075 where a distorted 
triangular (DT) lattice forms and the depinning is elastic. 



with the interstitial vortices depinning first. There is a 
strong enhancement of the anisotropy in the PP phase, 
with F-^/F^ reaching values as large as 5. For Fp < 0.175, 
the vortices form a distorted triangular (DT) lattice il- 
lustrated in Fig. Illf c) which becomes increasingly trian- 
gular with decreasing Fp. In the DT phase, only a small 
fraction of the pinning sites are occupied. Here the de- 
pinning transition is elastic and both the interstitial and 
pinned vortices depin simultaneously. The anisotropy of 
the depinning is lost and the critical depinning forces are 
isotropic in the DT phase. In general, F^ and de- 
crease with decreasing Fp] however, near the transition 
between the PP and DT phases, F^ increases with de- 



8 



creasing Fp as the depinning changes from plastic to elas- 
tic. We believe that this effect is similar to the peak in the 
critical depinning force observed for periodic [3a] and ran- 
dom arrays (39j| above the first matching field when the 
vortex-vortex interaction strength is varied. For periodic 
pinning, Ref. [s^ illustrated that the depinning force de- 
creases with increasing vortex- vortex interaction strength 
since the system depins elastically and all of the vortices 
begin to move at the same time. For weak vortex-vortex 
interactions, the interstitial vortices depin easily and flow 
plastically past the vortices at the pinning sites, so the 
depinning force decreases with decreasing vortex-vortex 
interaction strength in this regime. Between these two 
extremes, a peak in the depinning force occurs. A simi- 
lar effect is observed in vortex systems with very dilute 
random pinning arrays [39j . In Fig. llOf a) the vortex- 
vortex interaction strength is fixed; however, there is still 
a transition from elastic to plastic depinning. The peak 
appears only for x-direction driving, and this may be due 
to the existence of an easy shear mode in the a;-direction 
which is absent for y-direction driving. 

For B/B^ = 3, Fig. [TOjc) shows that pv and F^^ sat- 
urate at large Fp, while Fig. fTUT d) indicates that the 
anisotropy ratio also saturates at F^/F^ = 1.24. For 
Fp < 0.25, the depinning threshold decreases rapidly 
with decreasing Fp as the system enters the elastic de- 
pinning regime. We note that since the vortex lattice at 
B /B^ = 3 is triangular, there is no elastic energy cost for 
occupying the pinning sites and one-third of the vortices 
will always be located at the pinning sites for arbitrarily 
small Fp, unlike the situation at B/B^ = 2. There is a 
very small peak in near the plastic-clastic depinning 
transition. 

Figure [T2r a.b) shows the depinning thresholds and 
F^ as a function of Fp for B/B^ = 4 along with the 
anisotropy ratio F^ /F^. Both F^ and F^ saturate with 
increasing Fp, and the anisotropy for Fp > 0.25 is fixed 
at F^/F^ = 1.5. For Fp < 0.25, the depinning threshold 
decreases rapidly with decreasing Fp, and at the same 
time there is a drop in the anisotropy as the system 
passes from the plastic to the elastic depinning regime. 
In Fig. ^c,d) we illustrate F^ , Fv , and Fy/F^ for 
B/B^ = 5, where the vortex lattice lattice is disordered. 
Here the depinning is isotropic for all values of Fp. A 
peak in the depinning thresholds occurs for both direc- 
tions of driving near Fp = 0.25 at the transition between 
elastic and plastic depinning. 

We observe the same general trends for B/B^ > 5, 
including a saturation of the anisotropy with increasing 
Fp and a crossover from plastic to elastic depinning at 
low Fp. For higher values of Fp and Rp, multiple vortices 
are trapped at each pinning site. In this case, the vor- 
tices form vortex molecular crystals within the pinning 
sites [i^, which changes both the depinning threshold 
and the anisotropy. In general, the anisotropy is reduced 
when multiple vortex pinning is present. This means that 
the anisotropy is only observable in a regime where the 
pinning is strong enough to allow for plastic depinning 
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FIG. 12: (a) The depinning threshold in the a;-direction, 
(open squares), and y-direction, (filled circles) vs Fp for 
the triangular pinning lattice system in Fig. [5] at _B /B^ — 4. 
(b) The corresponding anisotropy ratio /F^ vs Fp. (c) F^ 
and F^ vs Fp for the same system at B/B^ — 5. (d) The 
corresponding anisotropy ratio FH /F^ vs Fp. Here the depin- 
ning is isotropic and a peak in the depinning thresholds occurs 
near the transition from the plastic to the elastic depinning 
regime. 



or channeling of vortices between the pinning sites, but 
not strong enough to permit multiple vortices to occupy 
each pinning site. 



B. Effect of Changing B^ 

We next examine the depinning thresholds and 
anisotropy at different matching fields for the triangular 
pinning lattice with fixed Fp = 0.85 but with increasing 
B^, achieved by increasing the density of pinning sites 
Up. The vortex- vortex interactions become more impor- 
tant for higher values of B^. For B/B^j, — 2, illustrated 
in Fig. [TWa.b). both F^ and FJ* show a peak feature, 
while the anisotropy remains near FJ'/F^ = 2.0. The 
depinning thresholds increase with increasing B^ for low 
B^ due to the fact that the increasing strength of the 
vortex-vortex interactions raises the interstitial pinning 
barriers. The threshold does not continue to monoton- 
ically increase with increasing B^ since the vortex con- 
figuration at B/Bcf, = 2 is a honeycomb lattice. As a 
result, the elastic energy cost of maintaining the hon- 
eycomb structure increases with increasing _B^, and for 
Bcf, > 0.2(f>o/X^, a portion of the vortices shift out of the 
pinning sites to form a distorted triangular lattice similar 
to that seen for low Fp in Fig. [TlT c). This causes a drop 
in the depinning thresholds with increasing i?^. We ex- 
pect similar behavior for other matching fields at which 
an ordered but non-triangular vortex lattice forms. 

In Fig. UM c) we plot F^ and F-^ versus B^ for the same 
system with B/B^ = 3, and in Fig. [THT d) we show the 
corresponding F^ /F'^. At this matching field, the vor- 
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FIG. 13: (a) The depinning threshold in the a;-direction, 

(open squares), and {/-direction, (filled circles) vs 
for a triangular pinning lattice system with Fp = 0.85 and 
Rp = 0.35A at B/B^ — 2. The vortex lattice has honey- 
comb ordering at low B^, but as B^ increases, the increasing 
strength of the vortex-vortex interactions causes some of the 
vortices to shift out of the pinning sites, (b) The correspond- 
ing anisotropy ratio F]/ /F^ vs B^. (c) F^ (open squares) and 
F^ (filled circles) vs B^ for the same system at B/B^ — 3. 
Here the vortex lattice is triangular, as seen in Fig.[S^d), so as 
B^ increases the vortex-vortex interactions become stronger 
and increase the value of the interstitial vortex depinning 
threshold, (d) The corresponding anisotropy ratio F^ / F^ vs 

Bj,. 
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FIG. 14: (a) The depinning thresholds F^ (open squares) 
and F^ (filled circles) vs for a triangular pinning lattice 
system with Fp = 0.85 and Rp = 0.35A at B/B^ = 4. (b) The 
corresponding F^ /F^ vs B^. (c) F^ (open squares) and F^ 
(filled circles) vs B^ for the same system at B/B^ — 5. (d) 
The corresponding F^ / F^ vs B^ showing that the depinning 
is isotropic. 



tex lattice is triangular as seen in Fig. [5]Jd). Therefore, 
the vortex positions do not shift as increases, unhke 
the case for B/B^ — 2, and the depinning thresholds in- 
crease monotonically with increasing B^. The anisotropy 
does not vary strongly with B^p. In Fig. [HTa.b) we show 
F^, and /F^ versus B^j, for the same system at 
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FIG. 15: The depinning thresholds F^ (open squares) and F]/ 
(filled circles) vs B^ for the system in Fig. [2] with Fp = 0.85. 
Lower curves: B/B^ = 0.5, showing that the depinning is 
isotropic and monotonically decreasing with increasing Fp. 
Upper curves: B/B^ = 1.0, where Fc — Fp. 



B/B^ = 4, where the vortex lattice is triangular as in- 
dicated in Fig. [ll^a) . The depinning thresholds increase 
monotonically with increasing B^ while the anisotropy 
remains constant at F}! /F^ « 1.54. At B/ B^ — 5, where 
Fig. md) shows that the vortex lattice is disordered, the 
depinning thresholds are isotropic and increase monoton- 
ically with increasing i?^, as illustrated in Fig. [HTcd). 
For any matching field B/B^, if B^ is increased above the 
range of values considered here, multiple vortex pinning 
by the pinning sites eventually occurs when the vortex 
lattice constant a becomes of the order of the pinning ra- 
dius Rp and very little distortion of the vortex lattice is 
required to shift the vortices into the pinning sites. The 
occurrence of multiple vortex pinning would alter both 
the depinning thresholds and the anisotropy. 



C. Anisotropy for B/B^ < 1 

For B/Bcf, > 1, interstitial vortices are present and the 
depinning threshold is determined by the strength of the 
vortex- vortex interactions, provided that the pinning is 
strong enough to produce a plastic depinning transition. 
For B/Bcf, < 1, the depinning threshold is controlled by 
a combination of the vortex-vortex interaction strength 
and the strength of the pinning sites. For a triangular 
pinning array at the matching and submatching fields 
of B/B^ — 1, 1/3, and 1/4, the vortex lattice is tri- 
angular and every vortex is trapped by a pinning site. 
In this case, the depinning thresholds are determined 
by the maximum force exerted by the pinning sites and 
are independent of the direction of the external drive, 
so the depinning is isotropic. This is in agreement with 
previous work, where an isotropic depinning threshold 
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FIG. 16: (a) The depinning thresholds (open squares) 
and (filled circles) vs T/Tm for a triangular pinning lat- 
tice system with Fp = 1.25 and B/B^ — 2. Here Tm is the 
temperature associated with the onset of vortex diflFusion at 
Fd = 0. (b) The corresponding anisotropy ratio F^/F^ vs 
T/Tm . (c) F^ (open squares) and F^ (filled circles) vs T/Tm 
for the same system at B/B^ — 4. (d) The corresponding 
FJ'/F- vs T/T^. 



was observed for vortices in triangular pinning arrays at 
B/B^ = l^^. In Fig. [m we plot F^^ and FJ' versus 
Bij, at B/B^ = 0.5 and 1.0 for a system with Fp — 0.85. 
The depinning thresholds are isotropic for both fields, 
and at B/B^ = 1, = Fp as expected. At B/B^ = 0.5, 
the vortex lattice is disordered since a triangular vor- 
tex lattice cannot match the pinning array at this fill- 
ing, as shown in Fig. 1 1 (c) of Ref . [41] . Although all of 
the vortices are pinned, due to the disorder of the vor- 
tex structure, some vortices experience stronger vortex- 
vortex repulsion from neighboring vortices than other 
vortices, and as a result the vortex-vortex interactions 
do not cancel at B/B^ = 0.5 as they do at B/B^ = 1. 
This lowers the depinning thresholds and causes both 
F'^ and to decrease with increasing due to the 
dependence of the depinning force on the vortex-vortex 
interaction strength. In general, we observe little or no 
anisotropy in the depinning thresholds for B/B^ < 1. If 
multiple vortex pinning occurs at B / B^ > 1, so that no 
interstitial vortices are present, we expect that the same 
type of depinning phenomena seen for B / B^ < 1 will 
appear and there will be little anisotropy. This suggests 
that the appearance of anisotropic depinning forces and 
the existence of voltage-current curves with different val- 
ues of dV/dFd in different directions are indicators of the 
presence of interstitial vortices in the system. 



D. Temperature Effects 

We next consider how robust the anisotropy is to ther- 
mal fluctuations. We anneal the system down to a finite 
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(b) 



FIG. 17: Vortex positions (black dots), pinning site locations 
(open circles), and vortex trajectories (black lines) just above 
depinning for a system with a square pinning array at Fp = 
0.86, Rp = 0.35A, and B^ = 0.0625(/>o/A^ at B/B4, = 2. (a) 
y-direction driving, (b) a;-direction driving. The same type 
of vortex motion occurs for both directions of drive. 
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FIG. 18: Velocity-force curves (T4) vs Fd for x-direction de- 
pinning (open squares) and {Vy) vs Fd for y-direction depin- 
ning (filled circles) in the square pinning lattice system from 
Fig. 1171 (a) B/B^ = 2, where the depinning is isotropic, (b) 
B/B^ = 4, with anisotropic depinning. (c) B/B^ = 12, with 
anisotropic depinning. 



temperature and then measure the depinning forces in 
the X and ?/-dircctions. The temperature is given in terms 
of the melting temperature Tm at which the vortices be- 
gin to diffuse significantly and the system is in a molten 
state. In Fig. [iga) we plot F^ and py versus T/T^ 
for a triangular pinning lattice system with Fp = 1.25 
at B/B^ — 2, and in Fig. [THT b) we show the corre- 
sponding anisotropy F^ /F^ . Both depinning thresholds 
monotonically decrease with increasing temperature and 
reach zero at T/Tm = 1-0. The anisotropy decreases 
with increasing temperature for low temperatures; how- 
ever, just below the melting temperature, F^/F^ passes 
through a peak. A similar trend is seen for B/B^ = 4, 
as shown in Fig. [TCT c.d). These results indicate that the 
anisotropy should be robust against temperature. This 
agrees with the experiments in Ref. 'si'l , which found that 
the anisotropy became more pronounced at higher tem- 
peratures. The enhancement of the anisotropy at higher 
temperatures could also result from the presence of weak 
random intrinsic pinning in the sample, which would re- 
duce the magnitude of the anisotropy but which can be 
washed out by thermal fluctuations. 
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FIG. 19: The vortex positions (black dots), pinning site loca- 
tions (open circles), and vortex trajectories (black lines) for 
the square pinning lattice system in Fig. [18] at B/B^ = 4. 
(a) The Fd = state where a triangular vortex lattice forms 
that is aligned in the y-direction. (b) The vortex trajectories 
just above depinning for driving in the y-direction, showing 
one- dimensional channels of vortices moving between adja- 
cent columns of pinning sites, (c) The vortex trajectories just 
above depinning for driving in the a;-direction. Here, every 
other column of vortices shifts in the y-direction so that every 
other vortex can join one-dimensional rows of vortices flow- 
ing between adjacent pinning site rows, while the remaining 
interstitial vortices are trapped between neighboring pinning 
sites, (d) The same as panel (c), but at a later time when the 
system has reached steady-state flow. 



IV. ANISOTROPY IN SQUARE PINNING 
ARRAYS 



For square pinning arrays, we find that the depin- 
ning thresholds at most fields are isotropic. This is a 
result of the fact that the same type of vortex motion oc- 
curs in both directions for most fields, as illustrated for 
B / — 2 in Fig. [171 In Fig.[l8Ua) we show the isotropic 
velocity-force curves at B / B^ = 2 for x and y-direction 
driving. Since the perpendicular directions of the square 
pinning array are identical, unlike the perpendicular di- 
rections of the triangular pinning array, it is not surpris- 
ing that most matching fields have the same depinning 
thresholds in both driving directions for the square ar- 
ray. Nevertheless, strongly anisotropic depinning occurs 
at B/Bff, = 4 and 12, as shown in Fig. [TSTb-c). 

In Fig. I19f a) we plot the vortex positions and pin- 
ning site locations for B/B^ — 4, where anisotropic de- 
pinning occurs. Here a triangular vortex lattice that is 



aligned with the y-axis forms. Under y-direction driving, 
the interstitial vortices fiow in one-dimensional channels 
between adjacent columns of pinning sites, as seen in 
Fig. [iW b) where 2/3 of the interstitial vortices have de- 
pinned. A simple channeling motion cannot occur for x 
direction driving, so is much higher than F}! . Just 
above the depinning transition for x-direction driving, 
as illustrated in Fig. I19r c). every other column of vor- 
tices shifts in the y direction in order to permit every 
other vortex to join a one-dimensional flowing channel 
passing between adjacent rows of pinning sites. The re- 
maining interstitial vortices remain trapped between the 
pinning sites. After this rearrangement, 2/3 of the inter- 
stitial vortices flow in the steady-state one-dimensional 
channels shown in Fig. [TQl d). which are similar to the 
channels that form for y-direction driving. The ground 
state shown in Fig. [TOT a) is degenerate, and a state with 
the vortex lattice aligned along the x-direction has equal 
energy. Thus, the realignment process seen in Fig. \Wi c) 
is simply a shift of the vortices into the other ground 
state prior to the onset of flow. Figure [TST b) shows that 
F'^ is approximately 1.75 times higher than F^ . Due to 
the symmetry of the square lattice, either the x or the 
y-direction can show a higher depinning threshold de- 
pending on the initial configuration of the vortex lattice. 
This is distinct from the triangular pinning lattice, where 
the higher depinning force always occurs in the same di- 
rection at a given field. If the annealing process for the 
square pinning lattice is repeated with different initial 
conditions, the vortex lattice has a 50% chance of align- 
ing with the x-direction, in which case the anisotropy will 
be reversed from that shown in Figs. [TH [TH If the driv- 
ing force is cycled, the velocity-force curve in the hard 
driving direction is hysteretic during the first cycle due 
to the realignment effect, while there is no hysteresis for 
the easy driving direction or for subsequent cycles in the 
initially hard driving direction. 

In Fig. [20] we illustrate the effect of changing Fp and 
B^ on the depinning anisotropy for the square pinning 
lattice sample with B/B^ — 4. We show a case where 
the initial vortex lattice orientation is along the y-axis, 
as in Fig. [TOTa). which gives Fj^ > F^. Figure ISOTb) in- 
dicates that the anisotropy saturates at F^ /F^ « 0.57 
for Fp > 0.075, while Fig. I^UT a) shows that the depin- 
ning thresholds also saturate above this value of Fp. For 
Fp < 0.075, the depinning is elastic and the anisotropy 
is lost. For fixed Fp = 0.85 the depinning thresholds 
increase monotonically with increasing B^, as seen in 
Fig. [20yc). while in Fig. [SOT d) the anisotropy passes 
through a shallow extremum of F^ /F^ = 0.5. 

We find a similar anisotropic depinning behavior for 
B/B^ = 12. The vortex configurations for this field at 
Fd = are illustrated in Fig. [2lT a). The vortex lattice 
is aligned in the y direction, but just as at B/B^ = 4, 
there are two degenerate ground states, and a vortex lat- 
tice that is aligned in the a;-direction would have the same 
energy. There are two columns of interstitial vortices be- 
tween adjacent columns of pinning sites. The overall vor- 
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FIG. 20: (a) (open squares) and (filled circles) vs 
for the system in Fig. [18] at B/B^ — 4 with the starting 
configuration shown in Fig. [T5] that produces a higher de- 
pinning threshold in the x-direction. (b) The corresponding 
anisotropy ratio F'^ifF^ vs F^. For small Fy the system de- 
pins elastically. (c) F^ (open squares) and F^ (filled circles) 
vs B^ for the same system at B/B^ = 4 and Fp — 0.85. (d) 
The corresponding F^ /F^ vs Fp. 



tex lattice structure is not triangular and there are dis- 
locations present in the lattice. Along certain columns, 
neighboring interstitial vortices in neighboring columns 
lie along a line tilted by —23° with respect to the x- 
axis, while along other columns, neighboring interstitial 
vortices lie along a line tilted by -1-23° with respect to 
the a;-axis. The dislocations in the vortex lattice are 
all aligned in the same direction, resulting in a smec- 
tic structure. This smectic state for the square pinning 
array has not been observed in previous work. There 
are two possible low-energy orientations for the vortex 
lattice, just as in the B/B^ = 4 case: the y-axis ori- 
entation shown in Fig. l^lTa). or the same state rotated 
by 90° and aligned with the x-axis. For depinning in 
the y-dircction. Fig. [HJb) shows that the two columns 
of interstitial vortices depin into flowing one-dimensional 
channels, while two interstitial vortices remain trapped 
between adjacent pairs of pinning sites. For driving along 
the x-direction, the same type of lattice reorientation 
found for B/ Bet, — 4 occurs at B/B^ = 12, as illus- 
trated in Fig. l^lTc.d). The columns of interstitial vor- 
tices shift in such a way that two one-dimensional rows 
of interstitial vortices form, while two interstitial vortices 
remain trapped between adjacent pairs of pinning sites. 
The final vortex configuration for x-direction driving, in 
Fig. EUd), is a rotated version of the configuration for 
y-direction driving seen in Fig. I2ir b). The same hys- 
teretic voltage-current response should occur for driving 
in the hard direction at B/B^ — 12 as at B/B^p — 4. 
In Fig. \Wic) the velocity-force curves at B/B^ — 12 
show an anisotropy /F^ k, 0.4. The more pronounced 
anisotropy at B/B^ — 12 compared to B/Bif, = 4 is 
due to the fact that the overall structure at B/B^ = 12 
is smectic and the dislocations are aligned in the y- 
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FIG. 21: The vortex positions (black dots), pinning site lo- 
cations (open circles) , and vortex trajectories (black lines) for 
the system in Fig.[T8]at B/B^ = 12. (a) The partially ordered 
vortex configuration at Fd = 0. The vortex structure has a 
smectic-type ordering in which all of the topological defects 
are aligned in the same direction, (b) The vortex trajecto- 
ries just above depinning for driving in the j/-direction show- 
ing two columns of interstitial vortices moving between each 
pair of pinning columns and two interstitial vortices trapped 
between adjacent pinning sites, (c) The vortex trajectories 
just above depinning for driving in the i-direction. A por- 
tion of the interstitial vortices shift in the y-direction in or- 
der to create ordered rows of flowing vortices and pairs of 
trapped interstitial vortices between adjacent pinning sites, 
(d) The steady-state channels that form for y-direction driv- 
ing. The general structural rearrangement of the interstitial 
vortices for depinning in the x-direction is similar to that seen 
at B/B4, = 4 in Fig. IMcdl. 



direction, further decreasing the depinning force along 
this direction. 

We expect that the depinning anisotropy in the square 
pinning arrays will be more difficult to observe experi- 
mentally than the anisotropy in the triangular pinning 
arrays due to the existence of twofold-degenerate ground 
states for the square system. If other forms of quenched 
disorder, such as intrinsic pinning, are present and the 
system is large, domains of the two different orientations 
could form which would render the depinning thresh- 
olds isotropic. On the other hand, we found that an 
applied drive can readily align the vortex lattice struc- 
tures in the driving direction at B/B^ = 4 and 12. 
Therefore, it may be possible to prepare the system in 
an aligned state using an external drive, and then mea- 
sure the anisotropy of the depinning forces starting from 
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FIG. 22: (a) (open squares) and (filled circles) for a 
honeycomb pinning lattice as a function of B/Bf , the honey- 
comb matching field, in a sample with _Fp — 0.85, -Rp = 0.35A, 
and Up = 0.194/A^ (b) The corresponding F^/F^ vs B/Bf 
shows that the anisotropy exhibits several reversals as a func- 
tion of field, (c) Fc (open squares) and F^ (filled circles) 
for a kagome pinning lattice as a function of B/B^, the 
kagome matching field, (d) The corresponding F^ /F^ vs 
B I B^ shows that several reversals of the anisotropy also oc- 
cur for the kagome pinning array. 



this aligned state. This procedure should permit the 
anisotropy in the square pinning lattice system to be ob- 
served experimentally. 



V. ANISOTROPY IN HONEYCOMB AND 
KAGOME PINNING ARRAYS 

We next examine anisotropy in honeycomb and kagome 
pinning arrays. Since the honeycomb and kagome arrays 
can be constructed by removing selected pinning sites 
from a triangular array, it might be expected that the 
anisotropy would follow the same trend found for the tri- 
angular pinning arrays. In particular, one could expect 
that the depinning threshold would always be higher in 
the y-direction. Instead, we find that the honeycomb 
and kagome pinning arrangements show an anisotropy 
that undergoes reversals as a function of the applied 
magnetic field. In Fig. [22{a) we plot and as a 
function of B/B^ for a honeycomb pinning system with 

Fp = 0.85, Rp = 0.35A, and Up = OAQA/X^. Here, 
B^ is the matching field for a honeycomb lattice con- 
structed from a triangular lattice with matching field 
B^, and we have B^ = 2/35^ [26]. Figure H^b) shows 
Fc / versus B^ for the honeycomb pinning array. For 
1 < B/B^ < 3, FV/F;^ < 1 and the depinning thresh- 
old is higher for x-direction driving. This anisotropy is 
reversed from that seen in the triangular pinning arrays. 
The reversal can be understood by examining the vortex 
positions at B^ — 1.5 in Fig. [53Ja). Each interstitial 



FIG. 23: Vortex positions (black dots), pinning site locations 
(open circles), and vortex trajectories (black lines) for the 
honeycomb pinning array in Fig. 1221 (a) At B/B^ = 1.5 
and Fd = there is one vortex per large interstitial site, (b) 
Vortex trajectories just above depinning for driving in the 
j/-direction at B/B^ = 1.5 showing one-dimensional paths 
between the pinning sites. At this field, F^/F^ < 1. (c) 
Vortex trajectories just above depinning for driving in the x- 
direction at B/Bf = 1.5. The interstitial vortices flow in 
winding paths around the occupied pinning sites, (d) The 
vortex configurations at B/B^ — 4 for Fd = 0. Here a 
portion of the interstitial vortices align in the x-direction. (e) 
Vortex trajectories just above depinning for driving in the x- 
direction at B/B^ — 4. In addition to moving channels of 
interstitial vortices, there are some interstitial vortices that 
remain pinned between adjacent occupied pinning sites, (f) 
Vortex trajectories just above depinning for driving in the y- 
direction at B/ Bf — 4 show a complex periodic pattern with 
all of the interstitial vortices moving. 



vortex is located at the position where the pinning site 
was removed from the triangular lattice in order to create 
the honeycomb lattice. There is a pin-free channel of mo- 
tion which the interstitial vortices can follow for driving 
in the y-direction, as seen in Fig. I^HT c). For i-direction 
driving, the path of the interstitial vortex is blocked by 
pinned vortices, creating a much stronger barrier for de- 
pinning, and once the vortices begin to move, they follow 
the winding paths illustrated in Fig. I^^Tb). 

In Fig. [M] we plot (T4) and (Vy) versus Fd for driving 
in the x and y-directions, respectively. The depinning 
in the honeycomb pinning array is anisotropic at both 
B/Bf = 1.5, shown in Fig. [24i;a), and at B/B^ = 4, 
shown in Fig.jMI^b). For 2.5 < B/B^ < 4 the anisotropy 
is reversed compared to the lower fields. This is due 
to the formation of n-mer states in the large interstitial 
spaces of the honeycomb lattice which permit a portion 
of the interstitial vortices to be aligned in the x-direction 
between the pinning rows. In Fig. [^5f d) we illustrate the 
Fd = vortex configuration at B/B^ ~ 4, where the in- 
terstitial vortices form triangular shapes within the large 
interstitial sites. The vortex at the top of each intersti- 
tial triangle forms a nearly one-dimensional channel in 
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FIG. 24: The velocity-force curves from the system in Fig. 1221 
Open squares: (Vx) vs Fd; filled circles: (Vy) vs Fd- (a) The 
honeycomb pinning array at B/B^ — 1.5 where F^/F^ < 1. 
(b) The honeycomb pinning array at B/Bf — 4 showing that 
the anisotropy has reversed and F]f /F^ > 1. (c) The kagome 
pinning array at B/B^ =4/3. (d) The kagome pinning array 
at B/B^ = 10/3. Here dVy/dFo < dV^/dFo, so there is a 
crossing of the velocity-force curves. 



the x-direction with the vortices at the base of the ad- 
jacent interstitial triangles, permitting easy depinning in 
the x-direction into the flow pattern shown in Fig. I^HTe). 
For depinning in the y-direction, where the threshold is 
higher, the complex but ordered flow pattern in Fig.[23Uf) 
appears. For = 4.5, the anisotropy reverses again, as 
seen in Fig. [22l 

Figure [22[c,d) shows that similar anisotropy reversals 
occur for depinning in kagome pinning arrays, where 
is the kagome matching field for a pinning lattice con- 
structed from a triangular lattice with matching field B^p 
and Bf = 3/4B^. For B/B^ < 5/3, Vig.^d) indicates 
that py /F^ < 1, but that this pattern reverses several 
times for increasing B/B^ . Just as in the honeycomb 
pinning array, in the kagome pinning array the reversals 
originate from the formation and alignment of intersti- 
tial vortex n-mer states in the large interstitial sites. In 
general, the anisotropy is smaller for the honeycomb and 
kagome pinning arrays than for the triangular or square 
pinning arrays. 

We illustrate the Fd ^ state for B/Bj^ = 4/3 in 
Fig. [25ja), where there is one vortex per large intersti- 
tial site. In Fig. EHb) we show the vortex trajectories 
just above depinning for driving in the y-direction at 
B/B^ =4/3. Each interstitial vortex diverts to the right 
or left around a pinned vortex, forming large asymmetric 
patterns of vortex flow around trios of occupied pins. The 
flow is slightly disordered, since in some cases two mov- 
ing interstitial vortices approach the same pinned vor- 
tex simultaneously, causing one of the vortices to move 
outside of the flow pattern temporarily. For x-direction 
driving at B/B¥ — 4/3, shown in Fig. ESf cl. the intersti- 
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FIG. 25: Vortex positions (black dots), pinning site locations 
(open circles), and vortex trajectories (black lines) for the 
kagome pinning array, (a) The Fn — state at B/B^ =4/3, 
where there is one vortex per large interstitial site, (b) The 
vortex trajectories just above depinning for driving in the 
y-direction at B/B^ =4/3 shows the formation of an asym- 
metric pattern that encircles groups of three pinning sites, 
(c) The vortex trajectories just above depinning for driving 
in the a;-direction at B/B^ = 4/3 shows a winding chan- 
nel of interstitial vortices between the pinning rows, (d) The 
Fd ^ state at B/B'^ = 10/3. (f) The vortex trajecto- 
ries just above depinning for driving in the a;-direction at 
B/B^ — 10/3 shows that a portion of the interstitial vor- 
tices move continuously through the system, (e) The vortex 
trajectories just above depinning for driving in the j/-direction 
at B/B^ = 10/3 shows a winding pattern of interstitial vor- 
tices. Here the vortex motion occurs in intermittent pulses 
rather than as a continuous flow. 



tial vortices flow around individual pinning sites in the 
sparse pinning rows, creating an asymmetric sinusoidal 
channel pattern. Figure IMTc) indicates that at this field, 
F^ I F^ < 1 and depinning is easiest along the y-direction. 
A similar anisotropy occurs for B/B^ = 10/3, as shown 
in Fig.[24i;d). In Fig.Iigd) we illustrate the Fd = state 
at B/B^ — 10/3, where the large interstitial sites cap- 
ture five vortices which form a pentagon structure. For 
depinning in the x-direction at this field, only a portion 
of the interstitial vortices move in winding channels be- 
tween the pinning rows while the large interstitial sites all 
capture three vortices, as seen in Fig. l^ST e). The intersti- 
tial vortices which are not part of the channeling flow still 
undergo a small circular motion as the flowing interstitial 
vortices move past. At higher drives, all of the intersti- 
tial vortices depin and a step appears in the velocity-force 
curve. For depinning in the y-direction at B/B^ = 10/3, 
shown in Fig. I25r f). the trajectories of the moving in- 
terstitial vortices are much more tortuous. Although a 
larger number of interstitial vortices spend at least part 
of the time flowing for depinning in the y-direction than 
for depinning in the cc-direction. Fig. I24f d) shows that 
dVy/dFD < dVx/dFD- This is because the vortex motion 
in Fig. HH^f) is not continuous. Instead, the interstitial 



15 



vortices hop by one lattice constant and then repin so 
that a pulse of motion passes though the system. For 
depinning in the x-direction, shown in Fig. [25je), the 
vortices are continuously moving through the system. 

VI. SUMMARY 

We have shown that when interstitial vortices are 
present for fields beyond the first matching field in reg- 
ular artificial pinning arrays, the transport response is 
anisotropic. For triangular pinning arrays, we find that 
the depinning thresholds are always higher in the y- 
direction, defined to be perpendicular to a symmetry axis 
of the pinning lattice. This is in general agreement with 
previous numerical and experimental studies up to the 
third matching field for triangular pinning arrays. We 
show that the anisotropy also occurs for higher matching 
fields and that at certain matching fields the depinning 
is isotropic. The degree of anisotropy can be controlled 
by changing the strength of the pinning. For weak pin- 
ning, when the depinning transition is elastic and all of 
the vortices depin simultaneously, the anisotropy is re- 
duced or destroyed. For strong pinning, the depinning is 
isotropic for matching fields at which the vortex lattice 
does not order. We find that the velocity- force curves 
can have different slopes for driving in the different di- 
rections depending on the number of interstitial vortices 
that initially depin. In some cases, although the depin- 
ning threshold is lower in one direction, the slope of the 
velocity-force curve is also lower in that direction, pro- 
ducing a crossing in the velocity-force curves for the two 
directions of driving. For fields at or below the first 
matching field, the anisotropy is strongly reduced, imply- 
ing that the anisotropy should disappear if multiple vor- 
tex pinning rather than interstitial vortex pinning occurs 



above the first matching field. The vortex dynamics can 
be distinctly different for the two driving directions, with 
ordered flow states occurring for driving in one direction 
and disordered flow states appearing for driving in the 
other direction. The anisotropy is robust against tem- 
perature fluctuations and can be enhanced near the vor- 
tex lattice melting transition. For square pinning arrays, 
the two perpendicular driving directions are symmetric; 
nevertheless, at certain matching fields where triangu- 
lar or smectic vortex structures form, a strong depinning 
anisotropy can occur. The easy-flow direction for square 
lattices is not fixed but depends on which of the two de- 
generate ground states is formed by the initial vortex lat- 
tice. A sufficiently large applied drive realigns the vortex 
lattice and sets the easy-flow direction in the direction 
of the drive. If the initial vortex lattice formed with do- 
mains of the two degenerate ground states, the domains 
can be eliminated and a pure ground state formed simply 
by sweeping the driving force. The anisotropic depinning 
behavior of this pure state can then be probed with small 
drives. Honeycomb and kagome pinning arrays have a 
smaller anisotropy than that seen for the triangular and 
square arrays. The anisotropy for the honeycomb and 
kagome arrays shows a series of reversals as a function of 
field due to the formation of vortex molecular crystals, 
which have orientational degrees of freedom that can lock 
to different angles. 
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